Cure to the Landau-Pomeranchuk and Associated 
Long- Wavelength Fermi-surface Instabilities on Lattices 
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Abstract 



The cure to the £ = 1 Landau-Pomeranchuk (L-P) instabihties in trans- 



(N 
> 
lO 

' lationally invariant fermions is shown to be a state with an anisotropic gap 



at the fermi-surface. For higher I and for fermions on a lattice, general cri- 
' teria for long wavelength instabilities and their cure are found in terms of 

' the derivatives of the single particle self-energy with respect to momentum 

for spin-symmetric instabilities and with respect to magnetic field for spin- 
antisymmetric instabilities. The results may be relevant to identifying hidden 
order parameters found in many metals. 



^ ' If the equilibrium distribution function n(k) of a translationally invariant (TI) Fermi- 



liquid is changed by 6n(k, a), 

Sn{k, a)=J2 [K' + K] Pe{e)S{e^ - /i), (1) 

e 

the system is stable only if 

F;''^> -(2£+l). (2) 

The superscripts s, a refer to the symmetric and antisymmetric spin-channels and F^'"' are 
the Landau parameters. These conditions derived by Pomeranchuk^are the generalizations 
of the positivity conditions on the compressibility, magnetic susceptibility and the specfic- 
heat, derived by Landau^. 
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Suppose the inequalities (3) are violated for some I. What is the symmetry and low- 
energy excitation spectrum in the state which cures the instability? Obviously, rotational 
symmetry must be lost for £ ^ 0. But rotational symmetry can be lost in at least two 
different kinds of states with quite different excitation spectra, states with a distorted Fermi- 
surface as is commonly presumed, or by a state with an anisotropic gap at the Fermi-surface. 
Starting from an equation of motion for the distribution function and using the relations 
between interactions and single-particle self-energy (Ward identitites) , I show that the loss 
of the rotational symmetry of the distribution function at least for the (. — 1 case must 
occur through a state with an anisotropic gap in the distribution function. For the spin- 
antisymmetric £ = instability, the answers from this approach are identical to that in 
the usual theory of itinerant ferromagnetism. For higher £ and for fermions on a lattice 
two classes of spin-symmetric instabilities are identified; one is cured by a state with an 
anisotropic gap at the fermi-surface, the other by a state with an anisotropic distortion of 
the fermi-surface. For spin-antisymmetric instabilities, the cure is a state with an anisotropic 
gap for spin-flip excitations at the fermi-surface. 

At the empirical end, there is evidence of gaps with "hidden order parameters''^"^ in 
many materials for which this work may be relevant. 

Let us specify quasi-particles created and annihilated by operators c^o-, Ckcr respectively 
with quasiparticle energy ek and an interaction function F(k, k', q, o", a') whose forward 
scattering limit q — > is the Landau interaction function. Specify further that F has a 
cut-off at so that it approaches zero for |ek — /^|, |ek' — A^l ~ ^c- Attractive Landau 
parameters might require that F(k, k', q, cr, a') come through the exchange of some Bosons 
with an upper energy cut-off generated by the repulsive bare interactions which have a 
much larger cut-off. (See discussion at the end.) 

Consider first the instabilities in the spin-symmetric channel. The equation of motion 
for small q is the homogeneous Landau transport equation: 

(ek+q/2 - ek-q/2 - U}) < C;^+q/2Ck-q/2 >s= (3) 
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[/(ek+q/2) - /(ek-q/2)] Ek' Fs{\i, k') < C+^q/2Ck'-q/2 >s, 

in which f{x) is the Fermi-distribution. 

If the forward scattering amphtude is attractive enough for any i to violate (2), no 
solution of Eq. (3) is possible for a; > in that angular momentum channel. The Fermi- 
liquid is then unstable and must reconstruct. The stability condition is given by putting 
cu — and then taking the limit q — > of Eq. (3) and may be written as 

$.(k+ q/2,k - q/2) > ±F,ik,k') ^^';'^''/'^ ' ^^'""-fh ^ik' + q/2,k- - q/2), (4) 

k' iek'+q/2 - ek'-q/2j 

[/(ck) - /(ekO]$.(k, k') = (ek - ek') < cj,Ck' >. . (5) 

The cure for this instability is to find a new set of quasi-particle operators c^^^, Ck,(T, with 
energies E^, which are linear combinations of the operators cj^, ^.Ck.o- such that the stability 

condition is satisfied. 

In the new (stable) state the stability condition is a new Landau transport equation 
analogous to Eq. (4), where again the hmit q — > is taken: 

l>,(k + q/2, k - q/2) > ± F(k, i,')lifK±S/A_^^^Lfl^^^i,' + q/2, k' - q/2), (6) 

k' l^k'+q/2 - £'k'-q/2j 

[f{Ey,) - f{Ey,^)]^,{k, k') = {Ey, - E^>) < cj^Ck' >. . (7) 

In Eq.(6), the sum is restricted to states such that -Bkiq are within about Vtc of the chemical 
potential jl. It has been assumed that the quasiparticle interactions in the new stable 
phase arc not different from those of the unstable phase. This is only true if the changes 
{E\f_ — ek) ^ Ef, i.e. in the weak-coupling limit. Strong-coupling corrections are not 
considered here. 

Replacing > by equality, the integral equation (6) serves as the self-consistency condi- 
tion to determine the temperature of the instability and the new quasi-particle energies £'k. 
Eq.(6)may also be derived in specific models and for fermions on a lattice. One such deriva- 
tion, which also serves to indicate a possible origin of the negative Landau parameters, is 
given at the end. 
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Suppose the instability occurs in the i-th angular momentum channel. Write F*(k, k') 
in a separable form and its i-th component as i/~^(0)F/P^(cos6'fc)-P£(cos^^), where z/(0) is 
the density of states at the chemical potential. Then the solution of Eq. (6) has the form: 

$,(k) oc PeicosOk), (8) 

and the the self-consistency condition in the limit q — > is 

1 - FMorY. panose,) ( ^^f>:iP ) • (9) 

Here , are the new one particle eigenvalues above and below the chemical potential 
respectively. 

The analogous equation when the instability is in the £-th spin-antisymmetric channel is 
1 = FMO)-E^/(cos^.) , (10) 

where E^pE^^^ are the new one particle eigenvalues for up and down-spin and above and 
below the chemical potential respectively. In the absence of spin-orbit scattering the spin 
may be quantized by the helicity, i.e. along and opposite k. 

Unlike BCS type theories which have an order parameter in an orthogonal channel, a self- 
energy in quadrature to the single-particle energy is not expected in the present problem, 
and 1 make the ansatz 

= ek + D^. (11) 

This ansatz is shown below to be consistent at least in the region just below the transition 
temperature. The form of D(k) is now picked through the microscopic theory for the approach 
to the instabihty. This is most transparently done for the i — and 1 channels. Consider 
the single-particle self-energy E(k, u) in the stable state as the L-P instability is approached. 
The instability in the symmetric £ = 1 channel is the vanishing of the effective mass m* 
given by^ 

m*/m = z-\l + {vfo)-'dE/dk\k„^)-' = 1 + Fi,/3. (12) 
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In Eq.(12) Vfo is the free particle velocity, and z is the quasiparticle renormalization ampli- 
tude, z = {1 — dY^/duj)~^. The first part of Eq. (12) is always vahd; the second part relies 
on TI. Since < 2; < 1, m*/m can be less than 1 only if E(k, ek) is negative/positive for 
k below/above k/, or on the energy shell below/above the chemical potential This is 
illustrated in Fig. (1) for the case the self-energy vanishes above Qc- Curve 1 is the unrenor- 
malized ek, curve 2 and 3 show the changes in the quasiparticle dispersion to ek + S(k, ek), 
as the instability is approached. At the instability, the dispersion at the chemical potential 
is vertical (i.e. m* = or divergent velocity) and negative beyond it. 

The cure to the L-P instability is to eliminate the region of negative dispersion in fig. 
(1) through an ansatz for the form of the self-energy, to be tested through solution of Eq. 
(9), which is a continuation of the form of the self-energy as the instability is approached. 
This choice is unambiguous, 

E> = + DCkf),E< = e^ - Dikf), for \E^- ii\< Dik^) 

^k = ek for \E^,,-ix\>n,. (13) 

In Eq. (13) D{kf) > all around the Fermi-surface, so that the ansatz, produces a (in 
general anisotropic) gap at the chemical potential. This ansatz used in Eq. (9) is sufficient 
to determine the transition temperature Tg and the magnitude and angle dependence of 
D{kf) just below it, as shown below. 

The alternative of deforming the Fermi-surface asks for a self-energy -D(k) which does 
not change sign at the chemical potential as one moves perpendicular to the Fermi-surface 
but changes sign on going around the Fermi-surface. This would require the ansatz 

D(k) = DvfSkpie), (14) 

where SkpiO) is the change in the Fermi-radius at 9, which one may seek to determine from 
the same self- consistency equation, viz (9). This ansatz is a continuation for the case that 
as the instability is approached, 

VpldE^uj, |q|)/5|q| |o,kp ^ -1, (15) 
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FIG. 1. Variation of the renormalized energy e(k) near the chemical potential in a particular 
direction due to increasing (1^3) attractive attractions with an energy cut-off 0,q. Curve labelled 
(4) is for Ff attractive beyond the L-P limit and therefore displays unstable negative dispersion at 
the Fermi-surface. 



6 



This corresponds to Fig — > oo in TI problems but even more generally to both Fqs, -Foa oo. 
The Mott or Wigner transition is a competing (first order) transition in this case. 

Consider next the £ = instabilities or the divergence of the compressibility or magnetic 
susceptibility: 

K m*/m X m*/m 



By employing Ward-identities^, one can show that 



(16) 



K 

- = (1 + -r^i^ + {vfor'dJ^/dk\,,)-\ (17) 

If the compressibility diverges without the effective mass diverging or independently of it, 
the L-P instability requires ^ oo- The ansatz for D to cure the instability is then again 
that it is discontinuous at the chemical potential just as in Eq. (13) leading to an isotropic 
gap in the single-particle excitation spectra in the stable state. In real physical systems, 
compressibility of the electronic system, ignoring its coupling to the lattice is generally not 
relevant. 

The (. — Q instability in the antisymmetric channel is of-course the prelude to ferromag- 
netism, which is signaled by Fq ^ — 1 or equivalently through Eq. (17) by cr • ^ oo. In 
this case the ansatz to be made in Eq. (10) is just the traditional ansatz in the mean-field 
theory of itinerant ferromagnetic instability, viz. 

£^k,a = Ck + (yD{kf) for |£;k,a - A*| ^ D{kf) 

Ek,a = ek for \Ek^a - A*| ^ f^c- (19) 

The adjustment of the up and down Fermi-surfaces so that the chemical potential for the 
up and down spins is the same is an additional condition which must be imposed here as in 
the traditional theory. 

For higher i, and for Fermions on a lattice, absent the Ward identities, only approximate 
derivation of the cure to the long wavelength instabilities are possible. (This is also the 



case of the familiar finite Q instabilities such as CDW or SDW.) One can rely on the fact 
that Eq. (3) can also be obtained in approximations such as the RPA or generahzed RPA, 
which use approximate interaction functions which we will continue to denote by F{k, k', q) 
and the approximate single-particle energies. The procedure to be followed is : Write the 
interaction in separable form in the irreducible representations (IR's) of the lattice. Check 
through solution of Eq. (3) or by some equivalent method if the interactions drives a mode 
in some IR unstable for some value of the coupling constant. Calculate the single-particle 
self energy near that point in a consistent approximation. On a lattice the self-energy may 
be expanded in irreducible representations of the lattice. ( I will continue to enumerate the 
IR's by i and denote their basis by 7^^ .) 

E(k,u;) = ^E^(|k|,a;)P^(cos^), (20) 

e 

For spin-symmetric instabilities, suppose that for any i, 

V dT,e{uj, \q\) / d\ci\\o,k^ ^ oo, or (21) 

^;^iaE,(a;,|q|)/a|q||o,kp^-l. (22) 

The single-particle energy is as exhibited in fig. (1) for the first case and has the opposite 
sign for the second case. The cure to the instabihty for the first case is found by the ansatz 
of the form of Eq. (13). The cure for the second case is found with the ansatz of the form 
of Eq. (14). The dependence of the functions -D(k) in either case is found from solving Eq. 
(10) as done below. An explicit example of this approach including the derivation of (21)for 
a specific case is given at the end. 

For the spin-antisymmetric instabilities, the RPA type equation show only one route to 
the cure. Suppose as the instability is approached, 

for some i. The cure is through an ansatz of the form (19). In the absence of spin-orbit 
scattering spin may be quantized by helicity. In that case or with spin-orbit scattering in 
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a crystal with inversion symmetry, no net moment is implied in any £ other than the £ — 
case by (19). The general solution in such cases has a gap function of the form Dij{kf)aj. 
The Magnetizaion direction rotates around the fermi-surface in the new state and there is 
an anisotropic gap single particle excitations with a given spin or for spin-flip excitations 

The relation of these instabilities on a lattice to the L-P instability except for the £ — Q 
or 1 cases is only conjectural. The instabilities are related to the divergence or zeroes of 
the eigenvalues of the generalized compressibilities d^F/duadup where Ua is the strain in 
the a-direction or the generalized magnetic susceptibihties d'^F/dHadHp, as are the L-P 
instabilities. 

I now return to the integral equation for D, Eq. (10) with the ansatz of Eq. (13) for 
arbitrary i for the TI or for the lattice case. Eq. (10) may be rewritten as, 

lp=i-Ffi/(0) / d9sm(9) 1,^, / deu(e) — -r- — -tvttt^- (24) 

^ ^' J ^' D{e) J-n, ^ 'l + cosh{(3e)sech{f3D) ^ ' 

The transition temperature to the new state is then calculated from Eq. (24) to be 

fir/2 

^ 1^ ( F^/{2t+l) \ ^^^^ 
l^l+F//(2f+l) J 

For T below Tg and {Tg -T)/Tg^l, Eq. (24) also gives that 

J desmeV]{cose)D\e) = 6{Tg - T)n^\F^\. (26) 

The requirement that D{9) > in Eq. (23) leads to 

D{e) ^ DoVficose), (27) 

where Dq may be determined from Eq. (26). (27) is not the only solution of the equation 
(26). Although all solutions with the ansatz (13) have an anisotropic gap at the chemical 
potential, the solution which minimizes the energy can be easily calculated to be the one in 
which the gap matches the kernel of the integral Eq. (26) most closely, as (27)does. 

Eq. (27) changes the one-particle density of states in a range of O(-Do) near the chemical 
potential /i. n must be determined anew to keep the number of particles conserved. For 
small Doi^(O), the new density of states for any £ is 
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viE) v{Q){E I Dof^ far E ^ Do (28) 

The specific heat in this state varies at low temperatures as T^^^ and the magnetic suscep- 
tibility as TV2. 

Similarly the spin-antisymmetric instabilities give rise to a gap for spin-flip excitations 
for any i with the value of the gap varying as (cos 9) . 

Landau or Pomeranchuk did not consider a cure for their instabilities. The symmetric 
spin-channel has been recently considered^' where the stable state is assumed to have a 
distortion of the Fermi-surface without a gap. A L-P instability for mesoscopic systems 
also has been proposed^. The solution similar to that here was proposed^ in connection 
with the special case of the pseudogap state of the cuprates. Several analytical-numerical 
calculations^°'^^ indicate a L-P instability in the "d-wave" channel in the Hubbard model. In 
analytic calculations collective spin-fluctuation modes have been argued to lead to negative 
Landau parameters'^ 

Derivation of Eq. (9) and of the form of the self-energy for a specific model of interac- 
tions: Consider bosons with bare energy fiq interacting with the fermions through coupling 
functions gf(k, k -|- q), 

Hboscms = E ^k&l^&k + ^(k, k + q)4+qCk(6q + ht^). (29) 
k k,q 

The nature of these Bosons is unspecified but it is necessary that /imq 0(7 (k, k + q) 7^ 0. 
Let the boson be unstable in the lim q — > due to the interaction with the fermions. The 
instability is cured by having a finite expectation values for 6q, 6+ in the limit q — > 0, which 
may be found from Eq. (9), 

^ g(k,k + q) + 

(^q) = - !L 7^ (Ck+qCk) (30) 

k ^"q 

The effective mean-field hamiltonian for the electrons is now 

Heff = J2 ^kCkCk + J2 ^(k, k + q)4+qCk, (31) 

k q 

V-(k, k + q) = ^(k, k + q)((6q) + (6lq)). (32) 
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ICk+qCk/ - ^; ^; wo; 



The one-particle Green's function for this hamiltonian in the self-consistent Brillouin-Wigner 
approximation is 

G'(k,a;) = ^— -, (33) 

where at a; = £'k = ek -l- -D(k, £'k), 

q ^k+q — ^k q 

The second equality in (31) follows from Eq.(21). Therefore, 

^(k, k + q)(/(£;k+qW(^k))((&q) + 

Inserting (6q), (61q) from Eq. (30), the self-consistency Eq. (9) is obtained with the Landau- 
function F(k, k') identified as lim q 0[— gf(k, k + c\)g{\s! , k' + q)/r2q]. Similar results can 
also be obtained from a purely electron-electron interaction model in which (6q -|- felq) in 
(26) is repalced by c,|^/+qCk' and c/(k, k -|- q) by a function V{k, k', q).In that case both the 
direct and the exchange interactions should be considered. 

If the Boson is a fluctuation of a non-conserved quantity, its propagator near the insta- 
bility is allowed the form, 

{iuj + a^q^ + e)-\ (36) 

e is the distance to the instability as a function, say of the couphng constants g at T=0. 
The fermion self-energy at the chemical potential for small deviation 5q from k^ normal 
to the fermi-surface is calculated to be oc \g{kFi ^F)P^Q'ln(q'c/(^9 -|- e) in 3 dimensions and 
oc \g(kF, kF)\'^{5Q + e)^^^ ln(gc/<^?) in 2 dimensions. (The singularities at k = k^T^ as a func- 
tion of energy are weaker.) In either case the derivative of the self-energy satisfies (21) 
as the instability is approached. The cure to the instability is therefore a state with an 
anisotropic gap at the chemical potential. This singularity in self-energy is reminiscent of 
the Hartree-Fock sigularity due to (unscreened) coulomb interactions. That singularities 
of-course disappears due to screening. Near an instability with a propagator of the form 
(36), the generated long-range interaction is protected. 
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